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Abstract 

The classical Lifshitz-Slyozov- Wagner theory of domain coarsening predicts asymp- 
totically self-similar behavior for the size distribution of a dilute system of particles 
that evolve by diffusional mass transfer with a common mean field. Here we consider 
the long-time behavior of measure- valued solutions for systems in which particle size 
is uniformly bounded, i.e., for initial measures of compact support. 

We prove that the long-time behavior of the size distribution depends sensitively 
on the initial distribution of the largest particles in the system. Convergence to the 
classically predicted smooth similarity solution is impossible if the initial distribution 
function is comparable to any finite power of distance to the end of the support. 
We give a necessary criterion for convergence to other self-similar solutions, and 
conditional stability theorems for some such solutions. For a dense set of initial 
data, convergence to any self-similar solution is impossible. 

Key words: Kinetics of phase transitions, domain coarsening, asymptotic behavior, 
self-similarity, stability, chaos. 
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1 Introduction 



The problem of Ostwald ripening is to explain the growth of the typical particle size 
observed in the late-stage aging of heterogeneously nucleated phase transitions. In the 
classical theory of Lifshitz and Slyozov® and Wagner, particles of the minority phase 
are assumed to be widely separated spheres that interact only through quasi-static dif- 
fusional exchange with a mean field. The chemical potential is proportional to curvature 
on the particle boundaries, and its gradient provides the flux that determines the particle 
growth rate. For background to the problem and a review of the physical literature, see 



refs. 15, 16 



In the dilute approximation with appropriate units, the volume v of any particle obeys 
the evolution law 

^- = K{v,e{t)) = v x 'H{t)-i, (i) 

where 9(t) is the same for all particles. The critical volume is 0(t)~ 3 = -yi? c (t) 3 , where 
R c (t) is the critical radius. Conservation of mass determines the value of R c (t) or 0(t), 
and if mass in the diffusion field can be neglected, total particle volume is conserved and 
one finds that the critical radius equals the average radius of currently existing particles. 
Particles with radius larger that R c {t) are growing, while small particles shrink to zero 
size in finite time and disappear. 

The LSW theory is concerned with a long-time regime in which most initially existing 
particles have disappeared but a large number remain. The system is characterized by 
the particle size distribution /(£, v), which we consider normalized so that f(t, u) du is 
the number of remaining particles with volume less than v, divided by the initial number 
at time 0. As particles of finite size are neither created nor destroyed, f(t, v) should 
satisfy the conservation law 

d t f + d v (A(v,e(t))f) = o, (2) 

where 

/•OO / poo 

0(t) = J o f(t,v)dv / J v 1 / 3 f(t,v)dv. (3) 

(These LSW equations also arise from a different model of particle clustering, the Becker- 
Doring equations, as recently shown by Penrose ©.) 

Lifshitz and Slyozov® and Wagner^llZP argued that for large times, the size distribu- 
tion should behave in a self-similar manner. They predicted: 

(a) The critical radius cubed grows linearly for large time, like 4t/9. 

(b) The particle size distribution approaches a form which is self-similar under scaling 
by the critical radius; and 

(c) In general, all scaled size distributions will approach the same self-similar solution. 
This solution is smooth, compactly supported and explicitly computable. 

The basic structure of the scaling predictions (a) and (b) have been validated in 
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experiments, but the form predicted in (c) was not found. (Also, the particular rate in 
(a) has proved difficult to measure.) Breakdown of the law (||) (which implies that all 
particles of a given size have exactly the same growth trajectory) is the most likely route 



for failure of the theory to explain experiment, as has long been recognized. See ref. 11 for 
a mathematical justification of this law by homogenization methods, that suggests that 
it is not necessarily valid, even in the limit of small volume fraction, if the electrostatic 
capacity of the particles does not also vanish. 

In this paper we address the mathematical problem of whether solutions of (||) in 
fact must have the asymptotically self-similar form predicted by the LSW theory. In the 
LSW analysis it is clear that the equations permit a family of self-similar solutions that 
are not infinitely smooth at their point of vanishing. However, the arguments of Lifshitz 
and Slyozov suggest that these are unstable. 

The present results indicate that the long-time behavior of the scaled particle size 
distribution depends sensitively on the behavior of the initial distribution near the end 
of its support. Let v(t) denote the supremum of the support of the volume distribution; 
we assume vo: = v(0) < oo throughout this paper. We prove, for example, that if the 
initial distribution has the property that 



i/> (u) := f(0,v)dv > au p (u > 0) (4) 

J VQ—U 

for some exponent p < oo and some a > 0, then convergence to the smooth self-similar 
solution predicted by the LSW theory is impossible. 

The non-smooth self-similar solutions are characterized by their exponent of vanish- 
ing. Given any exponent p in (0, oo) there is a unique (up to scaling) self-similar solution 
satisfying ipo(u) ~ au p as u — > 0. Our main result is to establish a condition on the initial 
data that is necessary for convergence to these self-similar distributions. The condition 
is a little weaker than requiring that lim u ^o ipo( u )/ uP exists — it states that the map ■i/'o 



must be "regularly varying at with exponent p" (see Definition |5.5| and the books of 
SenetaS 1 and Bingham et al.@). 

We conjecture that this necessary condition is also sufficient for convergence to a self- 
similar distribution characterized by a finite value of p. Numerical simulations that we 
have performed tend to support this conjecture. We prove it for initial data sufficiently 
close in a certain sense to some self-similar solution corresponding to a sufficiently small 
value of p. More generally, given arbitrary initial data that satisfy the necessary condition, 
we prove that the solution converges to self-similar form if the scaled inverse critical radius 
v(i) 1 / 3 6 l (t) converges as t — > oo. 

By making an arbitrarily small change to any given initial distribution in the vicinity 
of the tip of its support, one can easily arrange that the convergence criterion holds, or 
does not hold, at will. Our results then imply that there is a dense set of initial data 
for which convergence to any self-similar solution is impossible. If it is true that our 
necessary condition is also sufficient, then there is also a dense set of data that yield 
convergence to any given self-similar solution corresponding to a finite value of p. 

The basic idea underlying these results is a simple one. That is the notion that 
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the essential effect of time evolution is to stretch the size distribution, chop off the part 
corresponding to "dead" particles, and normalize to keep total volume constant. In other 
words, the action of the dynamics (up to a smooth distortion) is simply to "zoom in" 
as with a microscope on the details of the particle size distribution near the tip of its 
support. 

If the initial data behave in an irregular fashion near this tip on a microscopic scale, 
we can expect the solution to exhibit irregular behavior in time on a macroscopic scale, 
as particles are extinguished in order of increasing size and the distribution is normalized 
to conserve total volume. Indeed, the analysis in the proof of Theorem 5,10| below shows 



that, with respect to the logarithm of distance to the tip of support, the dynamics near the 
tip is rather well described by a continuous shift map. This suggests that the dynamics 
determined by the LSW model (Q) may actually be chaotic in a mathematical sense (with 
a dense trajectory and a dense set of time-periodic solutions). 

It would be wrong to infer, however, that the dynamics of Ostwald ripening is chaotic 
in any meaningful physical sense. The mechanism that produces sensitive dependence of 
long-time behavior on initial data in the LSW model depends heavily on the fact that the 
growth rate is the same for all particles of the same size. Lifshitz and Slyozov© already 
assumed some kind of breakdown in this rule in arguing that the non-smooth self-similar 
distributions would not be observed physically due to the effect of "encounters between 
different grains." It is plausible that a theory that seeks to quantitatively explain the 
self-similar distributions observed in Ostwald ripening should account for dispersion of 
growth rates among same-size particles. 

The notion that the asymptotic form of the size distribution in LSW theory depends 
on the precise mode of vanishing of the initial distribution has emerged recently (inde- 
pendently of the present work) in work of Meerson & Sasorov. Giron et al.,@ and 
Carr and Penrose.® The sole mathematical work in this connection is ref. ^, which con- 
cerns a simplified LSW theory for which the volume growth rate is modeled by v0(t) — 1, 
yielding linear characteristic equations. Carr and Penrose prove convergence to that self- 
similar solution characterized by ipo(u) ~ au p for initial data that have the property that 
Jo V'o^) du = A(x)x p+1 , where p > and A is differentiable with ^4(0) > 0. Further- 
more, they give examples of initial data for which no convergence to self-similar form 
occurs, and they also have results for data of infinite support. 

Giron et al.® expand a physical discussion given by Meerson & Sasorov$^P , and treat 
a class of models in which F(R,t), the number density of domains with radius R at time 
t, evolves according to a law of the form 

+ -S-(VF) = 0, V(R, t) = — ( — !— - - ) (n > -1). 
dt dR y ' ' v ' ; R n \R c {t) R) v " ' 

The critical radius R c (t) evolves so that the total volume J °° R 3 F(R,t) dR is con- 
served. These authors assert that a size distribution initially describable by a power 
law Ao(R m (0) — i?) A will for any t > have an expansion in which the leading term has 
the form A(t)(R m (t) — R) x , where the exponent A is invariant in time. It is then asserted 
that if a self-similar asymptotic regime is ever reached, it must correspond to the self- 
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similar solution having the same A 6 (0, oo). This assertion is supported by numerical 
computation and a linearized asymptotic analysis at the tip of support, for initial data 
having the form F (R) = £ x gi(0 + £ A2 #2(£) for £ = R m (0) - R > 0, where g\ and g 2 
are analytic and A2 > A > 0. (The condition A > corresponds to p > 1 in the present 
paper.) 

In the metallurgical literature, Brown (see refs. |^-|| and the references therein) had 
earlier asserted that an arbitrary initial distribution evolves to that self-similar solution 
"with the same tail," but by this he appears to always mean infinite-support self-similar 
solutions, for which the volume J °° vf(t, v) dv is infinite. Regarding distributions of com- 
pact support (which he calls "discontinuous" ) , he asserts they evolve to the distribution 
predicted by LSW. E.g., in ref. || he says "Both Hillert et al. (4) and Brown (5) have 
shown that an initial triangular distribution develops into the LSW distribution, albeit 
slowly." This kind of assertion is contradicted by the results we prove in this paper. 
However, Brown does make a number of persuasive arguments regarding the physical 
meaning of perturbing the distribution in the range of the largest particles in the system, 
and the ultimate dominance of the distribution of the largest particles. For example, in 
ref. |5] he mentions that the rate of coarsening can be slowed if the distribution favors the 
largest particles (as it does for self-similar solutions corresponding to small values of the 
exponent p). 

The structure of this paper is as follows. In Section 2 we reformulate the evolution 
law for the size distribution and recall the results proved in ref. 12 concerning the well 
posedness of the initial value problem. We show that the solution is given for finite time 
by an analytic distortion of the initial data. Section 3 contains results for long time that 
involve the unsealed variables. In particular, if the initial size distribution contains a 
Dirac delta at the tip of the support, the distribution converges to a stationary Dirac 
delta distribution as t — > 00, and if not, then v(t) increases to infinity. In Section 4, 
we classify self-similar solutions by exponent of vanishing. In Section 5 we scale the 
distribution function by the maximum particle volume v(t) and establish a variety of 
results as indicated above, concerning long-time convergence in the rescaled variables. 
We conclude in Section 6 with a discussion, and give examples of initial data yielding 
convergence and nonconvergence to self-similar form. 



2 Well posedness, regularity, stretching 
2.1 A well-posed formulation 

Some care should be taken in the physical interpretation of stability and instability results 
asserted for particle size distributions. One ought to consider what kinds of particle 
volume distributions are physically meaningful, and what kinds of perturbations make 
sense. It is plausible that in late-stage coarsening it is not possible to nucleate particles 
that are indefinitely large, for example. And it seems reasonable to imagine that small 
perturbations to the system would involve small changes to particle volumes. 

In this article we treat solutions for which the initial size distribution /(0, v) dv is an 
arbitrary probability measure with compact support. This admits distributions for which 
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a finite fraction of the particles all have the same size, for example. Compact support 
corresponds to the requirement that particle sizes initially are bounded. 

A theory for well posedness of the initial value problem for such solutions is developed 
in ref. |l2|. (For a family of closely related Lifshitz-Slyozov models, Collet and GoudonW 
have recently studied the initial-value problem for initial distributions /(0, •) that are 
integrable with finite first moment.) To discuss the results of ref. [l^, it is better to work 
with an equivalent formulation using different variables. First, let <p(t, v) denote the 
fraction of initially existing particles with volume > v at time t. With F(v) = 1 — ip(t, v), 
F is the usual distribution function for the particle volume distribution. When there 
is no danger of confusion we also refer to <p(t, •) as the "distribution function ip." The 
function v t— > <f(t,v) is decreasing and left-continuous at jumps with <p(t,0) = 1. (To 
be precise, we say (p is decreasing if <p{x\) < <p(x2) whenever x\ > X2, and similarly 
with increasing. A decreasing function need not be strictly decreasing.) The particle 
volume distribution is a measure, in general, given formally by f(t,v)dv = dF(v), that 
is, tp(t,x) = f(t,v)dv. 

Next we invert the relationship between v and ip and consider the volume v at time t 
as a function of ip £ [0, 1], taken so that <p v(t, ip) is decreasing and right-continuous at 
jumps with v(t, 1) = 0. For a finite system of particles with volumes ranked in decreasing 
order vo(t) > . . . > ujv-i(£), we have v(t, tp) = Vj for (p £ [jj, ^r). In general, 

v(t,x) = sup{y | ip(t,y) > x} for < x < 1 = max p. (5) 

We shall call the map (p \— * v(t,ip) the volume ranking for the system at time t, and 
refer to <p as the particle rank. The volume ranking determines the distribution function 
according to the prescription 



<p(t, y) = sup{x | v(t, x) > y} V 



( sup{x | v(t,x) > y] for < y < maxv(t,x), 
\0 for y > max-u(i, x). 

(6) 

The notion of distance between size distributions employed in the well-posedness 



theory of ref. |12| corresponds to the supremum-norm distance between the associated 
volume rankings. This is the least "maximal volume change" required to alter one size 
distribution into the other. Mathematically this corresponds to the notion of the L°° 
Wasserstein distance between probability measures, as shown in ref. 

The equation governing the evolution of the volume ranking v(t,ip) is just (|l|). In 



ref. 12, we proved the following theorem on the global well-posedness of the initial value 
problem. Let red be the set of right-continuous decreasing functions i>o: [0, 1] — ► R with 
■uo(l) = 0. On red we have the metric topology given by the supremum norm of the 
difference: 

\\v\ - v 2 \\ = sup \vi(ip) - v 2 (tp)\. 
f 

For any T > 0, C([0, T],rcd) is the metric space of continuous v. [0, T] — > red, with metric 
given by sup [0jT] \\vi(t, •) - v 2 (t, OH- 
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Theorem 2.1 Given Vq £ red, there is a unique continuous v: [0, oo) — > red and a unique 
9 £ L^ c (0,oo) such that 

/ v(t, <p)dp = / v ((p) dip 
Jo Jo 

for all t > 0, and 

v(t,<p)=vo(<p)+ [ (v(s,p) l ^9(s)-l)ds 



JO 

for all (t,(p) such that v(t,<p) > 0. Moreover, the map from vo to v is locally Lipschitz 
continuous from red into C([0, T], red), for any T > 0. 

The value of 9{t) is determined by conservation of volume to be 

/ r*® 

9(t) = p(t) v{t,y) 1/Z 'dip (7) 



o 

where <p(t) is the end of the support of v(t, •), i.e., 

(pit) = sup{0 G [0, 1] | v(t, 4>) > 0}. (8) 

The function <p is decreasing in time, but may have jumps. Consequently 9 can jump, 
but it is positive and bounded on finite intervals of time. For fixed tp, t i—* v(t,p) is 
Lipschitz continuous, and satisfies (|l]) for almost every t, as long as v(t,p) > 0. 

2.2 Stretching 

The result below indicates that for finite time the action of the dynamics stretches the 
volume ranking vertically in a smooth manner (by an analytic map). Thus, we can 
anticipate that whatever limiting behavior or irregularity is present in the initial data 
vq(p) in the limit ip — > will be preserved for any finite length of time. 



Proposition 2.2 Let (9,v) be a solution of as given by Theorem 2.1 . For x > 0, let 

V(t,x) be the solution of 

V{t,x)=x+ [ {V{s,x) 1/3 9(s)-l)ds 



Jo 

defined on the maximal time interval [0,t(x)) where V > 0. Then 

(a) V is analytic in x, i.e., whenever V(to,xo) > 0, the map x — > V(to,x) is analytic 
near Xq. 

(b) dV/dx is strictly increasing in time for each x. 

(c) v{t,p) = V{t,v {p)). 
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Proof: The formula in (c) follows by the uniqueness of solutions to the initial value 
problem for (|l|). 

Since the function w \— > w 1 ^ is smooth for w > 0, a standard proof shows that 
x i — ^ V(to,x) is continuously differentiable to any order, and each derivative is Lipschitz 
continuous in time. Almost everywhere in time, the derivative satisfies 

d_dV 8(t) dV dV 

dt dx (t ' x) ~ 3V(t, x) 2 /3 dx {tl xh dx ( ' x) ~ 5 



^-(t, x) = exp ( — ^ . ds, 



therefore 

dx \",-, -—*-J Q 3V ( s ,x) 2 /3 

and this is strictly increasing in time. This proves (b). 

To prove (a) we use a result related to a theorem of Bernstein. We recall that a 
smooth function g: (0, oo) — > M is called completely monotone if its derivatives satisfy 
(-l)V fe) > for all inte gers k > 0. As a consequence of Bernstein's theorem that every 
completely monotone function has a representation as the Laplace transform of a positive 
measure, it follows that every completely monotone function has an analytic extension 
to the open right half of the complex plane. (See ref. |l|, for example, for a proof.) 

We note that if V(to, xq) > 0, then it is easy to show that V(to, xq) < V(t , x) < oo for 
all x > xq. Thus, to prove (a) it suffices to show that if V(to, Xq) > 0, then the function 
x i ^ (dV/dx)(to, xq + x) is completely monotone. To do this, we show that the function 
given by 

u)(t, x) = —V(t, —x + xq) for < t < tQ, x < 0, 

satisfies d^u; > for all integers k > 1. This is true for k = 1 by (b). Suppose it holds 
for 1 < k < j for some j. Let G(y) = —(—y) 1 ^ 3 for y < 0, and note that the derivatives 
are positive for all k > 1. The function dito satisfies 

J^o, = 9(t)di(G(u)) = 9(t)(G'(u)dlu + R,(t, x)), 

where Rj(t, x) is a sum of terms of the form CaG^^^d^to . . . d" n uj which have positive 
coefficients and involve derivatives d^oj of order < j — 1. Then it follows that Rj(t, x) > 
and hence ^ ^ 

cV x iv(t,x) = J exp (J G'(u(r,x))dr \ 8(s)Rj(s,x) ds > 
for < t < to, x < 0. We obtain the desired result by induction. □ 



2.3 Evolution of the distribution function 

For the main part of our study of long-time behavior we will work with the distribution 
function (p(t,v), rescaled appropriately to study the development of self-similar form. 
To determine the time evolution of tp(t,v), we need to use (g) to get the appropriate 
generalized inverse function. The following result describes the evolution and implies 
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that it is a solution of the advection equation 

d t cp + A(v,0(t))d v ip = (9) 

whenever the initial data ipo = ip(0, •) is differ entiable. The characteristics of (^) are 
given by the function V of Proposition |2.2|. 



Lemma 2.3 Let ip(t,y) = sup{x | v(t,x) >?/}V0. Then 

f(t,V(t,y)) = ip (y) 

whenever < V(t, y) < v(t). 



Proof: Using Proposition 2/2 we find that 

<f{t,V(t,y)) = sup{x | V(t,v (x)) > V(t,y)} = sup{x \ v (x) > y} = <fo(y), 
provided < V(t,y) < v(t). □ 

3 Long-time behavior I 

In this section we derive some basic results as time tends to infinity without scaling the 
solution. Below we shall denote the maximum particle volume at time t by 

v(t) = v(t, 0) = supu (t, (p). 

Note that from (0) it follows 

v(t) 1/3 9(t) > 1, 

therefore v(t) is increasing. The total volume of particles in the system will be written 

V = / v (t, (p) dip. 
Jo 

3.1 The case of a Dirac delta at the tip 

First we consider the case in which a positive fraction of the particles have volume equal 
to the maximum possible in the system. In terms of the particle volume distribution / 
this means we consider measure- valued initial data /o that carry a Dirac delta at the tip 
of the support. That is, 

fo(v) = a5(v - v ) + f (v) 
where < a < 1, 5(v — vq) denotes the Dirac delta distribution at vq and 

poo 

lim / fo(y) dy = 0. 
v-^vo J v 
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For the volume ranking function v this assumption means that the solution v(t, ■) is 
constant on the interval [0,a). In this case we can characterize the asymptotic behavior 
as follows. 

Proposition 3.1 Assume that for some a £ (0, 1], the initial data satisfy vq(0) = vo((p) 
for < <p < a and vq(0) > vo((p) for <p > a. Then 

lim (p{t) = a, 

t— >oo 

and we have 

lim^) = (^ G ' °<^< a > 
t^oo K ^' \ 0, ip> a. 

lim ^ v) = l a > ^ v < V / a > 
t^oo^ y ' \ 0, v > V/a. 

Proof: From the assumption on the initial data it follows that for all t, v(t,ip) = v(t) 
for < (p < a and v(t, ip) < v(t) for (p > a. From conservation of volume it follows 



Correspondingly, 



poo 

V= v(t,(p) dip > av(t). 
Jo 



Hence v(t) is uniformly bounded, and since it is increasing the limit Vqo = lim^oo v(t) 
exists. 

Assume now that 

lim (p~(t) > a. 



t^oo 



(Since v(t) > V > 0, (p{t) > a for all t.) Then there exists ip > a such that v(t,ip) > 
for all t > 0. But then for almost every t > 0, using the inequality a — b> (a 3 — 6 3 )/3a 2 
for a > b > we have 

dt{v(t)-v{t,<p)) = 0(i)(t^) 1/3 -^) 1/3 ) 

eitwt) 1 ^ . .. . .. 

> ^?(»(*) -»(*.*>)). 

cincl hence 

v(t)-v(t,V>)>{v{0)-v (<p))^ 3V . 

Since by assumption v(0) — vo((p) > and since v(t) is bounded this contradicts the 
presumption that v(t,ip) > for all t. It follows that for each ip > a, v(t,<p) = for 
sufficiently large t, therefore lim^oo (p{t) = a. 

That = lim^oo v{t) = V/a follows by volume conservation. The limiting behavior 
asserted for <p(t,v) follows from the fact that a < <p(t,v) < (p(t) for < v < v(t) and 
<p(t, v) = for v > v(t). □ 
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3.2 The case of no Dirac delta at the tip 

Here and throughout the rest of this work we will assume that no positive fraction of the 
particles have the maximal particle volume v. For the initial volume distribution fo this 
means 



poo 

lim / fo(y)dy = 0. 



V— >1>o 

Proposition 3.2 Assume that vq(0) > vo(<p) for all <p € (0, 1]. Then 

(a) lim^oo v(t) = oo, i.e., the maximal particle volume tends to infinity. 

(b) liminfj^oo 9(t) = 0, meaning the critical radius is unbounded. 

(c) lim^oo (pit) = 0, i.e., the fraction of initially existing particles that continue to 
exist at time t decreases to zero. 

We do not know whether the result in (b) can be improved to say lim^oo 9{t) = 0. 
Proof: We prove that v(t) is unbounded using volume conservation. Suppose that 
v(t) < C for all t. As in the proof of Proposition 3A we find that for any tp > 0, as long 
as v(t, ip) > we have 

^)-<^)>(^(0)- Vo (^))e' /3C; - 

Consequently, for each ip, v(t, <p) = for sufficiently large t. This means that <p(t) — » as 
t — > oo. Then we infer that V < v(t)<p(t) — > 0, so V = 0, but the hypotheses guarantee 
V > 0. Hence v(t) is unbounded; since it is increasing, it tends to infinity. 

Suppose the claim in (b) fails. Then since 9(t) is bounded away from zero for t in 
any finite interval, there is some 5 > such that 

6{t) > 5 

for all t. From (a) and the fact that <p \— * v(t, ip) is right continuous, there exists e > 
and (to, <po) with <p$ > such that v(to, 'po) 1 ^ 3 ^ — 1 > e, which implies that 

^(t,p )=v 1 / 3 6(t)-l>e>0 

near t = to and so for almost every t > to- This yields that v(t,po) is unbounded, and 
this contradicts the inequality 

V= v(t,(p)d(p>(p v(t,(p ), (10) 
J o 

which is valid for all t since <p t— > v(t, (p) is decreasing. This proves (b). 
To prove (c), note that since v 1 / 3 < 1 + v we have 

<p(t) = o(t) [ v 1 / 3 (t,p)dp <e(t)(i + v). 



n 



By (b) there exists a sequence t n — ► oo such that 0(t n ) — > 0. Since <^(i) is decreasing the 
assertion (c) follows. □ 



4 Self-similar solutions 

For understanding the long-time dynamics, we shall study the distribution function or 
particle rank cp as a function of volume v according to the dynamical formulation de- 
scribed in Lemma ^H^. (Recall that the usual distribution function F(t, v) = 1 — (p(t, v).) 

We look for solutions in self-similar form. We begin with some general considerations, 
temporarily allowing for the possibility that the volume distribution has infinite support. 
Any self-similar distribution function ip for which total particle volume is finite and 
conserved must have the form 

cp(t,v) = a(t)-ip(a(t)v) (11) 
for some decreasing functions ip and a. The total volume satisfies 

V= (p(t,v)dv= / ip(s)ds, 
Jo Jo 

and we have 



Lemma 4.1 Any self-similar distribution function ip in the form (11) with finite total 
volume has compact support. 

Proof: Given the form in (0), by ( |I2| ) we may scale a and ^ so that #(£)a(£) -1 / 3 = 1. 
According to Lemma ^L^, being a solution requires that 

a(t)i>(a(t)V(t,y)) = a(0)$(a(Q)y) 

whenever < V(t, y) < v(t). For a nonconstant solution, a and ip should be Lipschitz, 
and differentiating we find that 

= a '(t)$ + s${s)) + a 2 {s 1/3 - l)^'(s). 

Separating variables, we find that a! jo? must be a negative constant. We set b = —a 2 /a', 
then a(t) = ba(0)/(b + a(0)t), and we must have 

= Vi(s) + (s + 6(l- S 1 /3))^( s ). (13 ) 
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As long as ip > we compute that 

_i (^ s A = r dx ^ 

\m) Jo x+ 6(1-^/3)- 

The function s *—> s + 6(1 — s 1 / 3 ) is convex with the value 6 > at s = 0, and at 
s = (6/3) 3//2 it achieves its minimum value 6(1 - ^46/27). This is positive if 6 < 27/4 
and nonpositive if 6 > 27/4. 

In the former case, we have if)(s) > for all s > and we compute that 

hm-lnfeUi±^Yl= / ( —r^- — \dx 



-0(0) V b J) Jo Kx + bil-x 1 / 3 ) x + b 



and this is finite and positive. Hence ip(s) ~ C/s as s ^ oo for some C > 0, and this 
means that the volume J °° i/j(s) ds diverges. 

In the case 6 > 27/4, on the other hand, the integral in ( |i~4| ) diverges as s — > so, where 
so is the smallest positive zero of s + 6(1 — s 1 / 3 ). (This zero is between 1 and the value 
(27/4) 3 / 2 /3\/3 = 27/8, and decreases to 1 as 6 — > do.) Then ip(s) — > as s — ► So, and we 
must put ?/)(s) = for s > So- D 

Thus we see that all self-similar solutions with finite total volume have compact 
support. To describe these solutions in a way that connects with the analysis to come, 
we rescale with respect to the tip of the support, introducing the variables u and ip*(u) 
that satisfy 

, V S i ( N V(t)ip(t,v) Sq4>(s) 

1 ~ U = W) = ~o' ^ (U) = ^^ = -^- (15) 
Note v = s/a(t) and v(t) = so/a(t) = v(Q) + tso/b. We set 

_ s + 6 _ 6s^ 3 
3so 3so 

then 6/so = 3k* — 1 and so = (3k*/ (3k* — l)) 3 , and k* increases from 1 to oo as 6 increases 
from 27/4 to oo. Furthermore we have 

s + 6(1-31/3) 

= u(k*Q(u) - 1), 

so 

where 

qm^ 1 -' 1 ;"' 173 ). (i 6 , 

The function Q satisfies Q(0) = 1, Q(l) = 3, and is given by the power series 

Q(u) = 1 + -u + —u 2 + . . . , 
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which converges for \u\ < 1 and whose coefficients are all positive. 
The self-similar profiles ip*(u) are increasing functions that satisfy 

T~ ln ^* = ( n) \ — TT ( 17 ) 
au u^KzQlu) — 1) 

for < u < 1, with the normalization 

/ 

Jo 

Note that for > 1, 

1 1 / re* \ — 1 



ip t (u)du = l. (18) 



u(k*Q(u) — 1) (re* — l)u \re* — 1 J u(n*Q(u) — 1) 



The last term is negative and remains bounded as u — > 0. Integrating fll7|), we find that 
self-similar solutions may be characterized according to the exponent p = l/(re* — 1), 
which governs the rate of vanishing as u — > 0. Henceforth, we write ^f p to denote the 
solution V* of (fiTh-dlD with re* = 1 + l/p- 



Lemma 4.2 For any p G (0, oo], equation ^) for the distribution function or particle 
rank has a self-similar solution determined by ( jl^ j and (15) with re* = 1 For 
p < oo £/ie profile has the form 

9 p (u) = a p (u)u p 

where a p is decreasing and analytic on [0, 1). Forp = oo (k* = I) we have ^oo(n) = o(u q ) 
as u — > /or a// g > 0. 

Equation (|l7]) can be explicitly integrated using the substitution 1 — u = x 3 and the 
expansion 

u(k*Q(u) — 1) = (a; — l)(x — d){x + a + 1) 



where 2a = — 1 + -y/3(4re* — l). For p < oo (re* > 1) we find 

(i - f ) P1 (i + 

where we have 



P2 



(19) 



3a 2 3(a + l) 2 



'' (o-l)(a + 2)' Pl (a-l)(2a + l)' P2 (a + 2)(2o + l)' 

In ref. |§| analogous expressions are given for the corresponding particle distribution 
$(r, x) := du*S(r, u) where 1 — u = x 3 . The solution singled out by Lifshitz, Slyozov 
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and Wagner is that corresponding to p = oo or = 1 and satisfies 



*oo(l) (l-x)5/3(l + a; /2)4/3- 1 > 



When interpreting these formulae, note that 

x = (I - n) 1 / 3 = 



^ R 



v(t) ) R(t) 

where R represents particle radius and R(t) is the maximal particle radius. 

As p —>■ in ( [19| ) we have a — > oo and we see that ^ p (u) — > 1 for < u < 1. That 
is, the particle distribution concentrates into a Dirac delta at the tip of the support. For 
the self-similar solution with profile the maximal particle volume satisfies 

v{t) = v(0) + Pl 



2p + 3 

Thus, as Brown® indicated, the rate of coarsening slows to zero as the particle distri- 
bution concentrates at the tip, which means p — > 0. 



5 Long-time behavior II 

5.1 Rescaled variables 

In the following we always consider the case that no positive fraction of particles has the 
maximal volume v. That is, we assume that the volume distribution carries no atom at 
the tip. To study the emergence of self-similar behavior for large times, it is appropriate 
to scale the volume by the maximal volume v(t). We introduce new variables via 

, v(t) , v , v(t)ip 

r = ln^, u = l-—-, <P=Arf L . 21 
v(0) v(t) V 

For each r > 0, the function u t— > tp(r, u) is left continuous and increasing on [0, 1] with 
V>(t, 0) = 0. With ipo(u) = v(0)ipo(v)/V , Lemma [2^ implies that the evolution of i()(r,u) 
is determined by the relation 

V?(t,«(t,u)) =e r VoN, (22) 
as long as <IA(r, u) < 1, where U(t,u) = 1 — V(t,u)/t'(t) satisfies 

— = (k(t)Q(U) - 1)U, W(0, u) = n. (23) 
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Here Q is the function defined in (^), and k(t) > | is determined from 

- -4 7T = 1 ^ = ,} \ I -(l-u)~ 2/3 ij(T,u)du. (24) 

The conservation of total volume is expressed by the identity 

/ 

Jo 



ip(r,u)du = l, (25) 



valid for all r > 0. 

For differentiable initial data, ip satisfies 



d T *P + (k(t)Q(u) - l)ud u ip = ifj (26) 



for < u < 1, r > 0. The characteristics of this PDE satisfy (p3|) . Note that k(t)Q(1) — 
1 > for all r, and therefore every point (ti,ui) £ (0, oo) x (0,1) lies on a unique 
characteristic that can be continued back to time r = 0. That is, iti =U{ti,u{) for some 
«l G (0,1). 

Self-similar solutions of @ correspond to stationary solutions of (12(f), for which we 
know k = > 1 is constant and V' = where p = 1/(k* — 1)- In general we are 
interested in investigating under what conditions ip{r,u) converges asr-> oo. We shall 
see that the long-time behavior of solutions is a sensitive function of the manner in which 
the initial data vanish as u —* 0. 

5.2 Tip behavior and the mean field 



Proposition 5.1 Let ip be determined by ( 2i )-(\23i), and suppose < p < oo. 



(a) If inf u>0 ipa{u)/vP > 0, then limsup^^ k(t) > 1 + |. 

(b) If sup u>0 ipo(u)/u p < oo, then h^mf r _>oo k(t) < 1 + ~. 



Proof: By Lemma 4.2, given the hypotheses of part (a), there exists ag > such that 



V'o(w) > ao^ p (u) for < u < 1. We aim to construct a "subsolution" of the form 

by taking a(0) = ao and ensuring that ip/ip- is increasing along the characteristics given 
by (|^) . By (^) it is equivalent to show the map r i— ► r — In (r, U (r, u) ) is increasing. 
Using (|l7]), we find that it suffices to ensure that a(-) is Lipschitz and for almost every 
r > satisfies 

<Aj) < 1 _ k(t)Q(u) - 1 = - k{t) 
a(r) ~~ K*Q(u) — 1 re* — 1/Q(u) 
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for < u < 1. Since 1 < Q(u) < 3 this is guaranteed if we require 

(k* - «(t))/(k* — 1), if - «(r) > 0, 
(«* - k(t))/(«* - §), if «*- k(t) < 0, (28) 
0, if ft* — k(t) = 0. 

With this choice, it follows tp(r, u) / a{r)^ p {u) > 1 for < u < 1 and for all r > 0. 

Suppose now that limsupK(r) < k*. Then liminf a' (r)/a(r) > and it follows 
a(r) — » oo as r — » oo. But this contradicts volume conservation. Hence limsupK(r) > /c*. 

Part (b) is proved in similar fashion. One constructs a supersolution of the form 
ip+(r,u) = a(T)^f p (u), ensuring that the inequality in (^7j) is reversed by interchanging 
the first two cases of (|2lf) . If liminf k(t) > then a(r) — > and again using volume 
conservation one obtains a contradiction. □ 



a(r) 



5.3 Necessary conditions for convergence 

Our aim in this section is to derive restrictions on initial data i/jq that must hold if 
convergence occurs in the rescaled variables. 

5.3.1 Consequences of convergence 

Lemma 5.2 Assume that lim T ^oo'/ , ( r ,'*) = 4'oo( u ) exists for each u E [0,1]. Then 
lim r _ +0O 1/k(t) = 1/Koo exists, where 1 < Kqq < oo, and ipoo(u) = \£ p (it), where p = 
1/(koo — 1) € [0, oo]. Moreover, as r — > oo, vJj(t,u) converges uniformly for u in any 
compact subset of (0, 1]. 

Proof: Since ip(r,u) is increasing in u, the dominated convergence theorem implies that 
lim r _ >oo 1/ K ( T ) exists. If the limit is zero, then since ipoo is increasing it follows from 
that Voo(n) = Voo(l) = 1 for all u G (0,1]. This is the case p = of the 

Lemma. 

If the limit is nonzero, then | < < oo. To show > 1, suppose otherwise. 
Then there is some uq > for which KooQ(u) — 1 < when < u < uq. For sufficiently 
large r this means that d T lA < if IA < uq; thus there is a characteristic satisfying 
< U(t, ui) < uq for t sufficiently large. Then (p2|) and the fact that ^(r, •) is increasing 
yield that ^>(r, it) > e T ^o(ui) for uo < u < 1, hence Jq 1 ip(r, u) du > e T (l— no)^o(^i) — > oo 
as t — ► oo. This contradicts volume conservation, so > 1. 

Volume conservation also implies that Voo(ito) > for some uo G (0,1). We aim to 
show that for any iti G (0, 1), i[ ? oa( u l) > and 

f n i (fit 

ln^oo(wi) - In ^00(^0) = T io := / —, tt- (29) 

First, note that if we define IAqq to be the solution of 

dU 

= (KooQ(Uoo) ~ l)Woo, ^oo(O) = U , 

or 
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then Woo(t) is denned for r in a maximal interval of the form (—00, To), with ^ QO (— 00) = 
0, U(Tq) = 1. Substituting n = U^t) into the integral in ( |29| ) reveals that Woo(no) = ni. 

Consider the characteristic for (^) passing through no at a large time To, as described 
by u = Uq{t) where 

= (k(t + r)Q(W ) - 1)W 0) ^o(0) = n . 

Then since k(t) — ► k m as r — > 00 and Q is smooth on (0, 1), standard arguments using 
Gronwall's inequality show that given any compact interval [a, b] C (—00, To), we have 

sup \14o(t) —^00(^)1 — ► as to — > 00. 

a<r<b 

Fix such an interval that contains and no in its interior. Note that it follows that 
d T Uo > for r 6 [a, 6] if to is sufficiently large. 

From (|22|) it follows that for any fixed r 6 [a, 6] we have 

•0(t o + r,W (r)) = e r ^(r ,u ) ->• e T V>oo(no) as r -> 00. 

If r < tio, then for sufficiently large To we will have Uq(t) < Z^oo(no) = ui, therefore 

t, ui) and it follows e T ip OD (u ) < ipoo(ui). Similarly, if r > no it 
follows e T ip 00 (uo) > ipoo{ui)- Therefore e T10 ^oo(no) = ipoo{ui), and this proves the claim 
in (|2£ 



It remains to show that ip^ is continuous at m = 1, so that (29) holds also with 
ui = 1. Since Q is increasing, for some S > we have 9 T Z^o > S whenever no < Uq < 1 
and To is sufficiently large. Let e > be small. Then since ip(r,-) is left continuous, for 
any sufficiently large To there exists no close to 1 such that iP(tq,uo) > V'oo(l) — e. With 
notation as above for this different value of uq we have Uq(—e) < 1 — Se, and therefore 

V>(t -e,l-5e)> V>(r - e,Uo(-e)) = e~ £ ip(T ,u ) > e~ e (V>oo(l) -e) 

The extreme members of these inequalities do not depend upon the choice of no, so we 
can take To — > 00 and deduce that Voo(l) < e e Voo(l — Se) + e. Since e is arbitrary, we 
infer that ip^ is continuous at 1, as required. 

From ([H]) and the normalizations in (25) and (18) it now follows that tpoo( u ) = ^p( u ) 
wherep= 1/(koo — 1), as claimed. Since ^ p is continuous and ip(r, ■) is increasing, asimple 
argument using pointwise convergence on a dense set in [0, 1] shows that ip(r, u) — ► ipoo(u) 
uniformly in u. □ 



As a consequence of Proposition [)7y and Lemma 5.2, we find the following result, 
which implies that ip(r, u) fails to converge for a large set of initial data. 

Corollary 5.3 Assume that lim r _ >00 i/j(t, 11) = ipoo(u) exists for each u £ [0,1]. Then 
there exists p G [0, 00] such that ipoo(u) = ^p( u ) an d 

inf ipo(u)/u q = for all q < p, 

u>0 
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and 

suprpo(u)/u q = oo for all q > p. 

These results imply the instability of the stationary solutions ^ p with respect to the 
sup norm topology, say, or any topology for which 

(a) the map from ip(r, •) to k(t) is continuous for fixed r, and 

(b) arbitrarily small perturbations of ^f p can yield initial data for which 

inf - > for some q < p. (30) 
u>0 ui 



Corollary 5.4 Given any p £ (0,oo], the stationary solution ^ p is unstable to pertur- 
bations that yield initial data that satisfy $3(\). In particular, the LSW solution (with 
p = oo) is unstable to all perturbations for which mf u> o^o(u)/u q > for some q < oo. 

See section 6 for examples and further discussion. 
5.3.2 Regularly varying functions 



Corollary |5.3| leads one to expect that if ip(r,u) — > ^f p (u) for < u < 1 and if p = 
— 1) < oo, then ipo(u) ~ cu p as u — > 0. This is almost but not quite exactly correct. 
In Theorem 5.10| below, we give a necessary criterion for convergence in this case that we 



conjecture is sharp. This criterion i s re lated to the concept of regularly varying function 
as is treated in the books of Senetafca' and Bingham et al.,® for example. 

Definition 5.5 A positive, measurable function g, defined on some interval of the form 
(0, a], is called regularly varying at (with exponent p G if 

lim f^f) = 1 for all A > 0. (31) 

X ^ Q + \Pg{x) y ' 

If p = 0, we say g is slowly varying at 0. 



Definition 5.6 We say that a real-valued, measurable function h, defined on some in- 
terval [A, oo), is locally linear at oo (with slant p £ M) if 

lim h(y + L) - h(y) = pL for all Lei. (32) 

y->oc 

If p = we say h is locally flat at oo. 

It is a classic result (cf. Lemma 1.1 of ref. [14| or Theorem 1.2.1 of ref. |2|) that if g 
is regularly varying at 0, then the limit in ([H|) is attained uniformly for all values of A 
in any compact interval [a, b] C (0,oo). Similarly, if h is locally linear at oo, then the 
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limit in ([32]) is achieved uniformly for all values of L in any fixed finite interval. This 
fact implies we have the following characterization of the two definitions just given. We 
define the oscillation of a function h on an interval [a, b] to be 

osc h(z) = sup h(z) — inf h(z) = sup h{z\) — h{zz)- 

z£{a,b] z£[a,b] z£[a,b] Zl ,z 2 e[a,6] 



Lemma 5.7 Suppose h(—lnx) = —lng(x). Then the following are equivalent, 
(i) g is regularly varying at with exponent p. 
(ii) h is locally linear at oo with slant p. 
(Hi) For all L > 0, 



osc (h(z) — pz) — > as y — > oo. 
ze[y,y+L] 



(iv) For all A > 1, 



1 as x —> 0. 



i*rf*e[*,Ax] 5 , ( z )/ zP 

The property of being regularly varying is invariant under a "nice" change of variables. 

Lemma 5.8 Suppose x t—* X(x) is C 1 near x = with X(0) = 0, X' > 0. Then, g is 
regularly varying at with exponent p if and only if go X has the same property. 

Proof: It suffices to prove the "only if" part, since we can write g = (g o X) o X^ 1 and 
X^ 1 has all the properties listed above for X. Assume that g is regularly varying at 
with exponent p, and let a < b with [a, b] C (0, oo). If A € (a, b) is fixed we have 

~X(X):=X X(XX) X - 



Ax X(x) 



as x — ► 0, hence A(x) G (a, 6) for sufficiently small x. Then, from the fact that (|3l|) holds 
uniformly for A G [a, b] we infer that 



g(X{Xx)) g(\{x)X{x)) ( \(x)Y 



1 as x — > + . 



A^(X(x)) A(x)p 5 (A») y A y 
This proves that go J is regularly varying at with exponent p. □ 



5.3.3 Convergence requires regularly varying data 

Recall from (||) that 

e T v{0) 
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where r = \n(v(t)/v(0)), and that V(t, •) is analytic in a neighborhood of u(t) with 
positive derivative. From Lemma and formula (^), we find the following. 



Lemma 5.9 For any r > andp E [0, oo), f/>(T, •) is regularly varying at with exponent 
p if and only if tpo has the same property. 

Here is our main result that gives a necessary condition for convergence in the rescaled 
variables to any one of the stationary solutions ^f p with < p < oo. (This result does 
not address the case in which the limit is the LSW-solution, the case p = oo in Lemma 

H) 

Theorem 5.10 Assume that for some p E [0, oo) we have 

lim ip(r, u) = ^p(u) for all u E [0, 11. 
Then is regularly varying at with exponent p. 

Proof: By Lemma [T^, to prove the theorem it suffices to show that for some tq > 0, 
V>( T o> ") is regularly varying at with exponent p. We consider first the case p > 0. Since 
k(t) — > as t — > oo, we may fix To > so that \k(t) — re*| < l/2p for all r > To- 
In the analysis to come we shall often make use of the variables 

y = — Inn, p = —lnip. (33) 
With po(y) = — hiV'o(^); w e have 

p(r,y(r J »))=p (y)-r (34) 
whenever y(r, y) > 0, where y = — InU satisfies 

^ = - K (r)Q(e~ y ) + l, y(0,y)=y. (35) 

For differentiable initial data, 

d T p-(K(T)Q(e- y )-l)d y p = -l. (36) 

In these variables stationary solutions are described via p*(y) = — In ^ (it) , and we 
have 

M( e -")-lK(l/) = l (37) 
for all y > 0. For < p = — 1) < oo it follows 

p*(y) = py - A *(y) ( 38 ) 

where lim^oo -A*(y) exists (and therefore A* is locally flat at oo). 
Define 

h(r,y) = p(r,y) - p*{y). (39) 
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The uniform convergence established in Lemma 5.2 implies that for any compact interval 
[a, b] C [0, do), lim T _ >00 h(r, y) = uniformly for y G [a, b]. 

To prove the theorem, it is enough to show that /i(to, •) is locally flat at oo. This will 
be done by following pairs of characteristics backward to time tq from large times f . The 
idea in the estimates that follow is that the total change in distance between such pairs 
can be bounded independently of f. 

Let y(r, y) denote the characteristic satisfying 



dy 

Or 



( K (r)Q(e~ y )-l), y(?,y)=y. (40) 



For f > t > tq we have 



hence 



dy 1 



y(r,y)>y + (f-r)/2p. (42) 



We also note that k(t) < Kq := + l/2p for r > to, so y(r, y) < y + 3Kq(t — r). 

Since Q'(0) = 1/3, there is some 5q > such that Q'(e~ v ) < 1 whenever y > 5q. 
Suppose that 5q < y < min(yi,y2)- Let 



y 1 (r) = y(r,y 1 ), y 1 =y 1 fo), 
Mr)=y(r,y 2 ), y 2 = Mto), 



(43) 



We compute that 

E(f,y 1 ,y 2 ):=(y 1 -y 2 )-(y 1 -y 2 ) = I K(r)(Q( e - y ^) - Q(e~ y ^)) dr. (44) 



TO 



Using d||) and the estimate ^(e"^ 1 ) - Q(e~ y2 )\ < e~ y , we find that 

\E{f,y u y 2 )\ < K e~y I e'^-^^dr < 2pK e~y . (45) 

J Tn 



From the relations (|3|), (|34|) and (g7|) we find that 

H(f,yi,y 2 ) := (h(r ,2/i) - h(T ,y 2 )) - {h(f,yi) - h(f,y 2 )) 

( ( Q(e- y ^)-Q(e- y ^) 
(k* — k(t)) = = dr. 46 

(«,Q(e-^M) - l)(K»Q(e-^«) - 1) 
As above we have the estimate 

\H\ < M (f) := p 2 e~y r |/c* - K(T)|e-( f_T )/* dr. (47) 

J T() 

Since k(t) — — > 0, it follows Mo(f) — > as f — > oo. 
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Let L > 0, and fix an arbitrary y > 5q. Put y 2 = y and choose y~\ > y~2 (depending 
on f) so that with y = y(f) = y 2 we have y + L = y\. By ( |45|) we have the bound 
Hi — V~2 < L + 2pKo, so we may deduce that h(f, y{) — > as f — > 00. 

Taking f — > 00 in (phf), we find that 

lim h(T ,y + L)-h(r ,y) = 0. (48) 

J/^OO 

Since L > is arbitrary, it follows that h(ro, •) is locally flat at 00. This finishes the 
proof of the theorem in the case p > 0. 

It remains to consider the case p = 0, when ^ p (u) = 1 for u G (0,1]. In this case 
p* = and in ( |46| ) we have H(f, 2/1,2/2) = 0. Moreover, hm T _ >0O «(r) = 00. We fix to 
sufficiently large so that k(t) > 2 for r > To- 

The strategy of the proof is the same as above. The key is to show that E(f, yij&te) 
is bounded independent of f for f > tq. Starting from (|40|), the first inequality in ( [41~1) 
remains valid, and moreover yields k(t) < —2dy/dr for r > tq. From (|4l) we get the 
bound 

\E(f,yi,ih)\< f K(r)e-WdT < 2 / e^ T )(-a^/ar)(ir < 2. 
Arguing as before, we obtain (f48|) and the claimed result follows. □ 



5.4 A sufficient condition for sufficiency 

The estimates in the preceding proof can be used to establish some other interesting 
results. We conjecture that the condition that tpQ is regularly varying at with exponent 
p > is actually sufficient as well as necessary for convergence to fy p as r — * 00. A 
general proof is lacking, but in the subsection to follow, sufficiency is proved when p > 
is sufficiently small and the initial data is close to in a sense related to the notion 
of regularly varying functions. The result given here reduces the general problem of 
sufficiency to proving the convergence of k(t). 

Theorem 5.11 Assume lim^oo k(t) = G (1>°°]« Then, linv-^cx) ^(^j u) exists for 
all u G [0, 1] if and only if tpo is regularly varying with exponent p = 1/(k* — 1). 

Proof: The "only if" part has already been proved. Suppose ipQ is regularly varying 
with exponent p. We first consider the case p > 0. Then for q < p we have ipo(u)/u q — * 
as u — ► and for q > p we have ipo(u)/u q — ► 00 as u — > 0. Since we assume k(t) — > 
as r — ► cxd, Proposition ^T] implies that p = 1/(k* — 1). For sufficiently large tq we 
have k(t) — 1 > l/2p whenever r > To, and we can use the estimates from the proof of 
Theorem |5. 10 . 

Our first step is to show that for any L > 0, 

sup \h(f, z) — b(f)\ — > as t — > 00, (49) 

z£[S ,6 +L] 
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where we define 6(f) = h(f, Sq). To prove this we use (plf), taking y 2 = Sq and y\ G 
[<5o, <5o + L]. Then J^i > 3^2 and it follows that E(?,y~i,y~2) < (the characteristics are 
diverging), so we have 

So + (f - T )/2p <V2<yi<V2 + L- 

Now @ yields that, 

|/i(f,yi) -6(f) | < osc h(T O ,z)+M (f)^0 

ze[y2,V2+L] 

as f — > oo uniformly for yi G [<5o ; + • This proves (|9|) . 
Next, we show that for any L > 5$, 

sup \h(f , z) — 6(f )| — > as f — » oo. (50) 

ze[o,L] 



To prove this, we use (|42|) taking y = 0, and use the identity ([5]) replacing t by 
n = f — 2p5o, which guarantees that y(ri,y) > 5q. Fix y~2 = Sq. For any y~\ G [0, L] we 
have 

So < Vi < L := L + 3£T (t - n) = L + 6p<5 ^o- 
Then since 1 < Q < 3 it follows from (^) that 

IM^ >2/i) — 6(f )| < osc_ h(ri, z) + 3p 2 / |k* — k(t)\cIt 

z£[5q,L] J n 



and this tends to zero as f —* oo uniformly for yi G [0, L], proving (|50|). 

Now ( |50| ) implies that for any L > 0, p(r, y) — 6(r) — ► /?*(y) as r — > oo uniformly for 
y G [0, L]. Since p(r,y) — 6(r) increases to oo as y — > oo for each r, we infer from the 
volume conservation identity (p^) that as r — > oo, 



/■OO /'OO 

e 6(r) = / e' p ^ y)+b{T) e- y dy ^ e~ p ^ e~ y dy = 1. 
Jo Jo 

Hence linv^oo 6(t) = 0. Then it quickly follows that lim^oo rf)(r, u) exists for each 
u G [0, 1]. This finishes the proof when p > 0. 

In the remaining case we assume lim,-^,^ k(t) = oo and ipQ is regularly varying with 
exponent p = 0. Then /i(to, •) = p(ro, •) is locally flat at oo for any tq > 0. We take To 
sufficiently large so that k(t) > 2 for r > to, and take ^2 = and y\ G [0, L]. Since 
(t) > 3M T ) > f — t for To < t < f , it follows E(f ,§1,1)2) < 0, hence f — To < y2 < 
V\ < 2/2 + Since -ff = in (f46|), we infer directly that 



sup |/t(f , 2) — h(f , 0)| < osc h(TQ,z)^0 

zG[0,L] z&[y2,V2+L] 

as f — > 00. With 6(t) = /i(t, 0), this means that for all L > 0, p{r,y) — 6(t) — > as 
t — > 00 uniformly for y G [0, L]. As above, from volume conservation we deduce that 
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6(r) — > and hence lim^oo ?/j(t, u) = 1 for all u £ (0, 1]. □ 
5.5 Stability and convergence for small p 

In trying to establish convergence of tp(r, u) as r — ► oo for regularly varying initial data, 
the main obstacle is achieving a priori control over k(t) — /c*. We can get such control 
when p is small, and obtain a stability result for the corresponding stationary solutions 
\I>p, with respect to perturbations that are small in a sense related to the notion of 
regularly varying functions. 

For perturbations that are small in the appropriate sense, the condition that the 
initial data is regularly varying is necessary and sufficient for convergence as r — > oo. 
In particular, the result below implies convergence to ty p (for p > small) for some 
initial data for which lim u ^o tpo(u)/u p does not exist. If the initial data are not regularly 
varying, convergence cannot occur. 

Recall that 

h(r,y) = p(r,y) - p*(y) = - ln(Y>(r, e ^)/* p ( e ^)). 
We define the flatness modulus of h to be 

vj(t, y) = sup qsc h(r,z). (51) 

y>y ze [5,5+1] 

Note that whenever < a < b < c, 

osc h(r,z)< osc h(r,z)+ osc h(r,z), 

z€[a,c] ze[a,b] ze[b,c] 

so for any positive integer n, 

osc h(r,z) < nw(r,y). (52) 

ze[y,y+n] 

Evidently, w(t, y) is decreasing in y, and we have w(t, y) — > as y — > oo if and only if 
h(r, ■) is locally flat at oo. 

Theorem 5.12 For sufficiently small p > 0, there exist positive constants 5* and Kq 
with the following property. If the flatness modulus of h satisfies w(0, 0) < <5*, i/ten 

tu(r,0) < K ro(0,0) (53) 

and 

sup |V(r, u) - %(u)\ < K tu(t, 0) (54) 

0<«<1 

for all t > 0. 

// m addition ipQ is regularly varying at with exponent p (equivalently, ifh(0,-) is 
locally flat at oo), then w(t, 0) -> ss t -> oo ; and linv^oo ip(r, u) = ^ p (u) for all 
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u € [0,1]. 

Proof: We start by getting a bound on \k(t) — k^\ where = 1 + l/p. Prom (^), using 
ip* = we have the representation 

y W.(i) ^(t,i); v ; 

= jT° (l - e h ^°)- fe (^)) e P*(o)-P*(y) (1 _ e -v)-2/3 e -v dy (55) 



«(r) 



We know that /0*(y) — /0*(0) > 0, and (|52; ) implies 

\h(T,0)-h( T ,y)\ < (y + l)w(T,0). (56) 
So, since |1 — e x | < Ixje^' for all x, the integrand in ( |55| ) is bounded by 

ro(r,0)(y + l)e ro(T ' 0)(y+1) (l - e^^/V 2 '. 
Thus there is a constant K independent of h and p such that if ro(r, 0) < | then 

1 1 



< i^ro(r, 0). 



(57) 



(58) 



k(t) k* 

It follows that if Kro(r,0) < 1/4^k 2 = p/4(p + l) 2 , then k(t) < 2k* and 

|k* - k(t)| < 2k 2 ^(t,0) < — . 

2p 

Next, we develop some a priori estimates for the flatness modulus, assuming it is 
small on some time interval. 

Lemma 5.13 For p > there is a positive decreasing function G p : [0, oo) — > R suc/i that 

POO 

/ G p (r)dr<Cp(p+l) 2 

JO 

for some constant C > independent of p, and the following holds. Assume that for 
some r* > we have Kw(f, 0) < p/4(p + l) 2 /or < f < r*. T/ien 

ro(f,0) <tp(0,f/2p)+ / G p (f - t)ot(t, 0) dr 



/or < f < t* 



Proof: We adopt the notation of the proof of Theorem |5.10| . First, let y > be arbitrary 
and suppose yi, J/2 £ + !]■ We use the identity (|46|) and the fact that 1 < Q < 3. 
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With any values of tq < f between and r* we have 

\h(f,yi) - h(f,y 2 )\ < |/i(r ,yi) - /i(r ,y 2 )| + 3p 2 / |k* - /c(r)| dr. 



Since the characteristics are diverging, yi and y 2 lie in [y, y + 1] where y = y(ro,y) > 
y + (f — To)/2p, and so since p 2 K 2 = (p + l) 2 , using (|58|) we find 



osc h(f,z)< osc h(ro, z) + 6(p + l) 2 K / ro(r, 0) dr 
[i/,s/+i] «e[M+i] Jro 



Taking the sup over y > y it follows that 



ru(f , y) < tu(r , y + (f - r )/2p) + 6(p + l) 2 K j w(t, 0) dr. (59) 



Next we use the same identity (Hq) but assume that y > 6q as chosen in the proof of 
Theorem 5.10| (so Q'(e~ y ) < 1 for y > 5q). Taking tq = 0, we obtain the estimate 



\h(f,m) - h{f,y 2 )\ < \h(0,vx) - h(0,y 2 )\+p 2 e-y / \k* - kIt)^-^^ dr. 

Jo 

From this we deduce that for y > So, 

w(f, y) < O7(0, y + f/2p) + 2(p + l) 2 e - 5 °i\~ / ro(r, 0) e - (f - r)/2p dr. (60) 



Now, we invoke ( [59] ) with y = and tq = max(0, f — 2p5o), and replace f in (pQ) by 
this ro and use y = 5q. Then, with 

rw-^ - J 6 Cp + i) 2 ^' 0<r<2p,5o, 
pl J " \ 2{p + l) 2 Ke-/ 2 P, r > 2pS , 

we find that the inequalities asserted in the Lemma are valid. □ 

Now we proceed with the proof of Theorem |5.12| . Assume p is sufficiently small so that 
Cp(p+l) 2 < i and assume Kzu(0,0) < p/8(p + l) 2 . Let M(f) = sup 0<T<f w(t, 0), then 
using Lemma p. 13 and a simple continuation argument, we find that M(f ) < 2ru(0, 0) for 
all f > 0. Thus we may take 5* = p/8K(p + l) 2 and K = 2 to prove the first estimate 
(|53|) asserted in the Theorem. 

The last part of the Theorem follows from Lemma |5.13j and the lemma below, together 
with the fact that if h(0, •) is locally flat at oo, then zu(0, y) — ► as y — > oo. One deduces 
that w(t,0) — ► as r — » oo. This implies /t(r) — > k*, and the convergence of ip(r,u) 



follows from Theorem 5.11 



Lemma 5.14 Suppose f, a, and g are positive bounded measurable functions on (0, oo) 
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and j Q °° g(s) ds = 7 < 1. A 



ssume 



f(t)<a(t)+ / g(t-s)f(s)ds 
Jo 

for all t >0, and a(t) — > as t — > 00. Then f(t) —* as i — > 00. 

Proof: Let F(t) = sup < s < t /(s), A = sup t > a(t). Then F(t) < A + 7 F(t) so F(i) < 
-4/(1 — 7) for all t > 0. Let F^ = limsup^^ /(£), and let e > 0. Choose to so 
/(s) < Foo + e whenever s > to- Then for t > to, 

f(t) < a(t) + f° g{t- s) ds + 7(^00 + e). 

1-7 Jo 

Taking the lim sup as t — ► 00, we find F^ < 7(Foo + e). Since e is arbitrary, it follows 
F oo = 0. □ 

It remains to prove the stability estimate (|54[). First, we get a bound on |/i(r, 0)|. 
Since Jq ip(r, u) du = ^f p (u) du = 1 we find 



1 



Mr,0)-h(r,y) 



1 ^(e-^e^ dy. 



Since ^ p (u) < a p (0)u p and |e x — 1| < Ixje' 21 ' it follows from (|56|) that provided ro(r, 0) < | 
we have 



^(t,0) 



1 



< 



a p (0)ro(r, 0) / (y + l) e ( y +1)m( - T ' 0) e~ (p+1)y dy < K x w(t, 0) 



where K\ is a constant. Provided we can guarantee vj{t, 0) < 1/2K\, we can ensure that 



e h(r,0) > 1 and then 



\h(r,0)\ = |lne h(T ' 0) -lne°| <2K 1 w(t,0). 



(61) 



From ( pq ) we then infer \h(r,y)\ < (y + i^) - ^^ 0) where F2 = 2K± + 1. Then it 
follows that for u = e~ y £ (0, 1], 



\^(r,u)-%(u)\ < 



-K T ,y) 



1 



V p (e- y ) < a p (0)\h(T,y)\e lh( - T ' y) \- py < K 3 uj(t,0) (62) 



where F3 is a constant, provided we can guarantee that tu(r, 0) < p/2. By (53), the 
bounds on ro(r, 0) that are needed to obtain ( |62"D will hold if 137(0, 0) < 5* for a sufficiently 
small value of 5*. This finishes the proof of the Theorem. □ 
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6 Discussion 

To recap and explain the main results of the last section, recall that u = 1 — v/v(t) 
and ip = (v(t)/V)ip where (p(t,x) = f°° f(t,v) dv, and suppose the initial data for tp is 
written in the form 

V>o(«) = (y) ¥X>fa>(l " «)) = ^He-" o(ln(1/u)) . (63) 

Here ipo must be increasing and is normalized so Jq ipo(u) du = 1. The results of section 
5 show that: 



(i) If the scaled distribution function ip converges pointwise as r — > oo, then the limit 
is ^ p for some p £ [0, oo]. 

(ii) If ipo{u) > au q > for some q < oo and a > 0, then ip cannot converge to $>oo, the 
solution favored by the LSW theory. 

(hi) If ip converges to ^ p for some p 6 [0, oo), then ipo is regularly varying at with 
exponent p. 

(iv) If p is sufficiently small and positive, and if the flatness modulus of ho is sufficiently 
small, then ip converges to ^> p if and only if ipo is regularly varying at with 
exponent p. 

Recall that i/jq is regularly varying at with exponent p if and only if the function 
/io in ( |63[ ) is locally flat at oo. can be replaced by u p in ( |63[ ) and this statement 
still holds.) Thus, if ho is bounded (implying ^o(^) > au p ) and is not locally flat at oo, 
then ip cannot converge pointwise to any limit as r — > oo. For a concrete example, any 
periodic or quasiperiodic function Hq will suffice, e.g., ho = a — bs'my, corresponding to 
ipo(u) = Cu p e bsinXn ^/ u \ 

If ho is C , the mean value theorem implies that the flatness modulus of ho satisfies 

ro (y):=sup osc h (z) < sup \h' (y)\, (64) 

y>y z£[y,y+l] y> y 

and so ho is locally flat at oo if h' (y) — > as y — > oo. As an example, suppose ho(y) = 
— ln(l + ey) + c, so ipo(u) = C^ p (u)(l + e\n(l/u)). Then ho is locally flat at oo and 
^o(O) < £. So if p is small and e is small, the corresponding solution ip converges to ^ p 
in this case, despite the fact that ipo(u)/u p — > oo as u — > 0. 

Sharp representation formulae are known for regularly varying functions. For exam- 



ple, from Lemma 1.3 of ref. |l4j, it follows immediately that ho is locally flat at oo if and 



only if there exists b > such that for y > b, 



ho(y) = h(y)+ f h{t)dt, 
Jb 
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where f\ and /2 are bounded measurable functions such that lim^oo fi(y) exists and 
limj / ^ 00 / 2 (y) = 0. 

The results we have established in this paper show that the asymptotic behavior 
of solutions is exquisitely sensitive to the behavior of the initial data near the tip of 
the support. The conditions in (ii) and (hi) that ensure nonconvergence in the rescaled 
variables depend only on the initial size distribution of an arbitrarily small fraction of 
the largest particles. That is, they depend only on the initial volume ranking vo(ip) for 
values of the particle rank ip in an arbitrarily small interval (0,£o). 

It is clear that one can always perturb the sizes of an arbitrarily small fraction of the 
largest particles, by a maximum amount as small as you like, to guarantee that ipo is not 
regularly varying at and ipo(u) > avP for some p. This means that, with respect to the 
topology given by the sup norm distance between volume rankings (a natural topology 
for well-posedness of the initial value problem), there is a dense set of initial data that 
yield nonconvergence to self-similar form. 

By similarly perturbing the sizes of the largest particles, one can guarantee that i^q 
is regularly varying at with any given exponent p. If it is true, as we conjecture, that 
having such data implies that V converges to the self-similar profile *$> p , then there is a 
dense set of initial data that yield convergence to ^f p , for any p £ [0, oo). 

One can also make an arbitrarily small perturbation to make the new volume ranking 
constant on a small interval of the form (0, e). This means that the size distribution 



contains a Dirac delta at the tip, so by Proposition 3.1 the perturbed size distribution 
will converge to a stationary Dirac delta as t — > oo. 

We have not discussed in this paper any conditions on initial data for convergence 



to the smooth self-similar profile favored by the LSW theory. Lemma |5.2| implies 
that k(t) — ► 1 is necessary for such convergence. In this case, characteristics near the 
tip spread more slowly and different sorts of estimates seem to be needed to study the 
convergence question. 
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